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The t-J z model is the strongly anisotropic limit of the t-J model which captures some general 
properties of the doped antiferromagnets (AF) . The absence of spin fluctuations simplifies the an- 
alytical treatment of hole motion in an AF background and allows us to calculate the single- and 
two-hole spectra with high accuracy using regular diagram technique combined with real-space ap- 
proach. At the same time, numerical studies of this model via exact diagonalization (ED) on small 
clusters show negligible finite size effects for a number of quantities, thus allowing a direct compar- 
ison between analytical and numerical results. Both approaches demonstrate that the holes have 
tendency to pair in the p- and d-wave channels at realistic values of tj J. The interactions leading 
to pairing and effects selecting p and d waves are thoroughly investigated. The role of transverse 
spin fluctuations is considered using perturbation theory. Based on the results of the present study, 
I . we discuss the pairing problem in the realistic t-J-like model. Possible implications for preformed 

' pairs formation and phase separation are drawn. 
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I. INTRODUCTION 



a 
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Q ' The physics of holes moving in an AF background has received much attention because of its possible connection 
to high-T c superconductivity. A microscopic realization of this physics is given by the t-J model which was first 
I ' introduced as a conjecture and then derived as a low-energy limit of a realistic model representing the electronic 
structure of cuprates ||J|]. Recent comparison of the t-J model results to the angle-resolved photoemission data |4j 
■^J- ■ showed that the overall shape of the experimental single-hole band can be fitted satisfactorily using the pure t-J 
model. To account for the detail line shapes, one needs to include more distant hopping terms (f , t", etc.) J|||. The 
presence of these terms also follows from a careful mapping study Q. Numerical studies of the t-J model by means of 
ED and other methods suggest the presence of hole binding in the physical parameter range tj J ~ 3 [0-[9|, and that 
the dominant symmetry of the pairing correlation is d x 2_ y 2. This finding also supports the common belief that some 
variant of the t-J model is able to describe the physics of the real compounds. However, to show that this d-wave 
binding in the t-J model is relevant to the physics of high-T c materials, one has to clarify how strong and generic the 
■ reasons behind this pairing are. Unfortunately, numerical studies alone are of limited use for this purpose, so that an 
\ analytical study is necessary to develop an insight into the problem. 
I ■ In this paper we attempt such a study by investigating the Ising limit of the t-J model - the t-J z model. It is 
well-known to be a simplified limiting case of the t-J model, but only recently has its properties been fairly well 
understood. In the pioneering work by Bulaevskii et al (!(]] the single-hole energy in the tj J <C 1 limit has been 
calculated. Subsequently, Brinkmann and Rice |pd| considered the J = limit and introduced the retraceable-path 
approximation. More recent works fli"2|-]l6| used the self-consistent Born approximation (SCBA) to study the t-J z , 
t-J, and t-t'-J models. Using this approximation, an analytical expression for the hole Green's function in the t-J z 
model for arbitrary values of tj J has been found by Starykh and Reiter ]T7| ). The pairing of holes in the t-J z model 
has been considered in the works by Trugman [|3~8|] and by Shraiman and Siggia p9[ , where important results were 
obtained. Some other studies used different modifications of the variational approach ]2C| , |2l| ] and obtained similar 
results. However, as we shall show, previous analytical treatments involved approximation which prevented one from 
obtaining the correct answers. 

The t-J z model has also been studied using numerical methods |t]]. Barnes et al p^] compared their ED results 
on a 16-site square lattice with some analytical results, and Riera and Dagotto [^3j developed a modified Lanczos 
technique which enable them to study the one- and two-hole ground states of the model on lattices as large as 50-sites. 
(For a review see Ref. 0). It was shown that finite-size effects in the ED data of the t-J z model are much smaller 
than those of the t-J model p3| . 

In this paper we develop an analytical treatment of the t-J z model based on the results of Ref. Comparison 
with our ED data allows us to justify the validity of this analytical approach. Furthermore, we use the results of this 
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study to shed light on the problem of whether the pairing in the ci-channel is a generic feature of t- J-like models, and 
what interaction defines the symmetry of the ground state of the system. 

Our numerical ED results for the t-J z model are obtained on a 32-site square lattice using the same method as in 
Refs. This is up to now the only ED results on the largest lattice. This numerically exact method allows 

us to analyze various properties of the single-hole ground state as well as the two-hole bound states with different 
symmetries. 

Our paper is organized as follows. Sections II and III describe the representation and Feynman rules used in this 
paper. In Sec. IV we solve the single-hole problem and compare the results to ED data. In Sec. V the two-hole 
problem is considered in detail. In Sec. VI we summarize our results and draw conclusions. 



II. REPRESENTATIONS OF THE HAMILTONIAN AND OPERATORS 



The Hamiltonian of the t-J z model is 

H = - l E (4**r + H.c.) + JJ2 [SIS* - \NiNj] (1) 

(%j)o 

where the summation runs over nearest-neighbor bonds (ij), $ ia — c\ a (l — n i& ), ci a = (c| CT )^, N t = cj^c^ + cj^Cjj. 

Since we are modeling electrons on a 2D square lattice, the coordination number z = 4 and the spin S = 5. 

Assuming the presence of long-range antifcrromagnetic order one can introduce a representation for the Hubbard 
operators (operators of spin and constrained fcrmion) through the spinless fermion and boson operators. The choice 
of the representation is, in general, arbitrary and is motivated by the problem to be solved. In our work we use a 
generalization of the Dyson-Maleev (DM) representation for the Hubbard operators. As usual, such a representation 
conserves the algebra of the operators but extends the Hilbert space of the problem. 

In sublattice A = {| } the DM representation is given by 

St = 1(XP X?) = \- a\ ai -\fth> Sr= = a\(l aja, S\ = X? = a t 

2,T = *f =//(!- ajoi) , 5l T =*F = fi, *l = = f \ a . _ ^ = x io = f . a i (2) 

Ni = 4^1 + ifri = X P + X P = 1 - x i° = 1 - flfi ■ 

In sublattice B = {J,}, the corresponding representation is given by 

% = l( X P - X P) = -\ + b fa + ^j<?i , Sj=X?=b j , S ]=X^ = b](l-b]bj- g] 9j ) , 
c j l=X?=g]b j , c^=X]°= 9j b], ~ Cji =Xf = 5 t(l-b) & .) , c) i =Xf=g j , (3) 

4t^t + = X P + X P = i- x !° = i- 9)9, ■ 

Thus, we have two types of fermions and bosons associated with the A and B sublattices. The Hamiltonian in Eq.(jl]) 
expressed in terms of the new variables Eqs.(^,||) is 

TL => H.q + TL\ + H.2 , 

n, \j Y. + % + flfi + g}gj] , 

(i,3) 

n 1 ^tY^[fl9j(b]+a i ) + R.c.] , (4) 

(ij) 

Ui [fUiM* + sUiatfyi] - \ J E H h h + "»< + n H + 2 «MM • 

(ij) (i,j) 

Here we have divided the Hamiltonian Eq.(|4|) into three parts: TIq, consisting of the linear spin-excitation and hole 
terms, Tli consisting of the bare hole-magnon interaction term, and TI2, consisting of the nonlinear hole-magnon, 
static hole-magnon, direct hole-hole, and magnon-magnon interaction terms. The third and fourth terms in TLq take 
into account the energy of the four AF bonds broken by a hole introduced into the system. In the following this 
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energy is included in the ground state energy and the corresponding terms are omitted. Note that we have omitted 
a constant term. The same Hamiltonian Eq.(^|) in k-space is given by 



q k 

Hi = 7k- q [/L q <?k^ + 9Lj^l + B.C.] , 

k,q 

H 2 =-4t 7k- q [/il_ q _ qi+q2 .9k& q a qi a q2 +^_ q _ qi+q2 / k a q 6 qi & ' 

k,q,qi,q 2 

- 2J 7q-qi[/k- q + qi /k& q & qi +gLq+ qi .9kaq a qi] 

k,q,qi 

" 2J 7q/k-q5k'+q5k'/k-4J ^ 7qOqi -q 6 q 2 +q 6 ' 



Q2j 



(5) 
(6) 



- , q 2 Ll qi ) 



k,k',q 



where 7k = (cos(fc x ) + cos(fcj,))/2, and all summations are restricted to the first magnetic Brillouin zone (BZ). 



III. PROPAGATORS AND VERTEX FUNCTIONS 



The Feynman rules for the model in Eq. (g) are: 

(i) Bole propagator: G°(e) = G° f{g) (e) = [e + iO}- 1 

(ii) Magnon propagator: D°(w) = D°, b) (uj) = [w - 2 J + iO}- 1 

(iii) Bolc-magnon vertex: 

r°(fc, q,k- q) = 4i7k- q ■ £ei,£-w 



te 



a,b 



a(b) 



fig) 



k,e 




k,£ 



q,co 



g(f) , fig) 
k-q, el k-q, el 




The indices indicate that / (or g)-holes can emit only a (or 6)-magnons and absorb only b (or a)-magnons. This is 
an important feature of the hole-magnon interaction, and is a consequence of spin conservation. In the following we 
omit the indices. We use the following shorthand notations — 8^^-^ denotes 8{e\, e — uj) and ^ w means |H. 
(iv) Bole-two-magnon vertex: 

qLcol q,co 



F°(k,q,qi,k - q + q x ) = -2J7 q _ qi ■ S eue - U+U)1 



a(b) — 




fig) 
k-q+ql, el 



Note that / (or <7)-holes interact with a (or £>)-magnons only, 
(v) Bole-three-magnon vertex: 



r%(k,q,qi,q2, k-q-qi + qi) = -4£7 k _ q -5 ei!£ _ w _ Wl+W2 



a(b) 
fig) 



q2,co2 



q,co 



b(a) 




ql.col 

- a(b) 

- g(fl 



k,e k-q-ql+q2, el 



Note that there is no vertex diagram for the reversed order of emission- absorption because of the non-hermitian nature 
of the corresponding terms in the Bamiltonian Eq. (ra) . 
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(vi) Hole-hole vertex: 




where we use the notations k = (k, e) and q = (q, u>) for the momentum and frequency of the hole and magnon 
respectively. 

Thus the model Eq.(||) involves fermion-fermion (hole- hole), boson-boson (magnon-magnon), and three types of 
fermion-boson (hole-magnon) interactions. In general the problem of finding the low-energy excitations in such a 
model is very complicated. However, as we will show below, in the case of the t-J z model it is possible to take into 
account the result of renormalization almost exactly using a regular diagrammatic treatment. 



IV. SINGLE HOLE 



If we ignore 7^2) the Hamiltonian in Eq. (g) for a system with a single-hole is similar to the problem of a single 
electron interacting with the local phonon mode. When the interaction (t in our case) is small compared to the energy 
of the boson (2J), we can use perturbation theory and consider only the lowest contribution to the hole energy, which 
corresponds to the virtual emission-absorption of a single magnon. At larger t many-magnon intermediate states 
must be considered. Such a problem, if solved self-consistently, must take into account the renormalization of the 
fermion-boson interaction, or crossing diagrams. However, in the t-J model it is well known that the lowest correction 
to the bare hole-magnon vertex T° is exactly zero since the hole-magnon interaction conserves the spin, or more 
exactly, the pseudospin associated with the sublattice type. Therefore, the crossing diagrams can be omitted and the 
renormalization of the hole energy is given only by the series of diagrams where the magnons are absorbed in exactly 
the reversed order in which they arc emitted. This is the so-called self-consistent Born approximation (SCBA), which 
in the case of the t-J z model is identical to the "retraceable path approximation". This coincidence is due to the local 
nature of the spin-excitations in the model. The corresponding diagrammatic equation for the hole Green's function 
is shown in Fig. 111. Note that the vertices are unrenormalized. 




FIG. 1. SCBA for the hole self-energy. Bold lines represent the dressed Green's functions. 

The retraceable path approximation implies that the hole is confined by a "string" of spin excitations and the hole 
motion is completely incoherent, i.e. the hole form a localized state. 

The first non-zero contribution to the hole self-energy from the crossing diagrams is shown in Fig. 0. It is of sixth 
order in t and is small compared to the contribution of the "retraceable" diagrams of the same (sixth) order due to 
the geometrical factor. 
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FIG. 2. Trugman's diagram. The numbers represent the sites where a magnon is emitted (absorbed). 



This diagram is the first member of the family of the so-called "Trugman paths" . Their contribution to the hole 
energy is small in a wide range of t/J, but they are responsible for the finite coherent band of the hole. Real-space 
consideration shows that such a diagram corresponds to the hole's motion around an elementary square loop in which 
the hole makes one and a half turn. As a result all spin excitations are cured and the hole is translated along the 
diagonal of the square plaquette, i.e. to a next-nearest neighbor. As we have noted, the corresponding correction to 
the energy is small and we will neglect these diagrams in the rest of this Section. 

Since G (e) and D°(uj) are momentum-independent and the hole- magnon interaction depends only on k — q, the 
full Green's function G is also k- independent. Thus the equation shown in Fig. [j] can be simplified as 

E(c) =^G(e- W )D°( W ) E( r k-q) 2 = ~ 2 J ) > 

uj q 

G(e) = [e-4t 2 G(e-2J)}- 1 . (7) 
The solution of Eq. (Q) at all t/J has been found in Ref. j[7). One can try an ansatz for G of the form: 

This transforms Eq. (Q) into a difference equation: 

Y(e + 2J)+Y(e-2J) = ~^Y(e) , (9) 
which is the recursion relation for the Bessel functions. Thus G(e) is given by 

where J u (x) are the Bessel functions. 

Such a form is the consequence of the continued fraction form of the Green's function, which in turn is the result 
of the retraceable path approximation. Note that in Ref. jl7j J z — 2 J. Thus the poles of the Green's function in the 
SCBA approximation are defined by the zeros of the Bessel function. Since all these zeros are real, the hole spectral 
function consists of the set of delta-function peaks corresponding to the energy levels of the quasiparticle states in 
the "string" potential well. Inclusion of the Trugman processes does not change this result significantly because the 
resulting bandwidth is much smaller than the separation of the levels. 

Although the solution @ of Eq. (g) is an almost exact solution of Tio + 7ii, it is not a solution of the original t-J z 
Hamiltonian. Besides the crossing diagrams, there are other contributions to the self-energy originating from Ti.2- In 
order to compare with other analytical approaches and ED, it is necessary to account for these contributions. 

To make this statement evident let us now consider the two limiting cases t/J< 1 and t/J 3> 1, and compare the 
results of the SCBA with the known facts from the real-space approach and ED numerical data. 

For t/J -C 1 the ground-state energy for the hole is given by Eqs. (Q) as 

e SCBA = _^ (n) 

This is inconsistent with the result of the real-space approach 

eo = -—. (12) 
3J v ' 
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The reason for this discrepancy is that the spin-flip (magnon) has lower energy (3 J/2 compared to 2 J) if it is located 
in the neighborhood of the hole, and the Hamiltonian Ho + Hi does not take this difference into account. Note that 
throughout this paper we assign zero energy level to the state with one static hole, i.e. Eq = E I a smg + 2 J = 0. 
According to Ref. fnj, in the large t limit the ground-state energy given by Eq. ( |To| ) is 

e scBA = _ 4t + 2 /3 t(J/t) 2 ^ - 2 J , (13) 

where /3q = 2.34 is the first zero of the Airy function. There is no alternative asymptotic expression in this limit which 
is accurate up to order J, but the first two terms which are of order t and t 1 ^ 3 J 2 / 3 can be compared to the results of 
the "string" approach by Shraiman and Siggia (SS) jl9| and the fitting of 16-site ED data by Barnes et al [^2] 

e$ s = -2V& + 2.74t( J/i) 2/3 , (14) 
e§ D = -3.63i + 2.93i(J/<) 2/3 . 

This comparison makes one to suspect that SCBA (Eq.(|l3|)) overestimates both the absolute value of the depth of 
the "band" and the slope of the walls of the linear confining potential. 

The origin of these discrepancies becomes evident when we consider the real-space picture of the retraceable hole 
movement. Fig. ^| shows the real-space images of three components of the spin-polaron wave function: the one-, two-, 
and three-magnon strings. The first important observation is that the coordination number in the creation of the first 
magnon is indeed z = 4, but is z — 1 = 3 in the creation of each next magnon. This is why the factor in the first 
term of the SS expression Eq. ( |l4| ) for the hole energy is 2^3 (instead of 4). If t is much larger than J, the average 
length of the strings contributing to the spin-polaron wave function is large. Then the weight of the components with 
coordination number different from others becomes insignificant. Since this component is the "bare" hole, it also 
means that the quasiparticle residue of the state decreases as t grows and in the t -> oo limit the spectrum will be 
incoherent. This result was first obtained by Brinkmann and Rice flil[ in the J = limit. 
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FIG. 3. One-, two-, and three-magnon strings 

Having in mind that the energy of spin excitation is determined by the number of broken AF bonds associated with 
its creation, one can see that the energy of the first magnon in the string is 3J/2 (3 bonds, J/2 each), but the energy 
of each subsequent magnon is just J (2 bonds). This is true for all strings of length / = 2 and for most other longer 
strings except for those which have self-tangencies. Since the latter cases are relatively rare due to the geometrical 
factor, one can assign the same energy I ■ J + J/2 to every string of length I. A more accurate account for the energy 
of the self-tangient strings makes only tiny changes to the energy of the whole excitation. Thus, making the changes 
2 J — > J and z — > z — 1 in the expression for the SCBA hole energy will change the factor 2/?o = 4.68 in the second 
term of Eq. (0) to i 1/& (3 Q = 2 .81, leaving it in reasonable agreement with Eq. (|l4|). 

Summarizing the results of the comparison one can say that the reason for the failure of the SCBA is that Ho +Hi 
is not the same as Ht-j and that an improved formalism must take into account the following facts: (i) the energy 
of the first magnon in a string is 3 J/2, (ii) the energy of the nth magnon (n > 1) is J, (iii) the coordination number 
is z — 1 = 3 for I > 1 and z — 4 for I — 0, where I is the length of the string. After some thinking using the above 
real-space arguments one can suggest the following modification to the equation for the Green's function (Q) to meet 
all these requirements 

G(e) = [e-4t 2 G a (e-3J/2)}- 1 , 

G a (e) = [e-3< 2 G Q (e- J)]" 1 , (15) 

G ( e ) = 1 J -E/42%/3t/J) 

V3t J_ e/J _ 1 (2x/3t/J) ' 

Such a modification gives the correct hole energy in the small t limit. In the t — > oo limit the resulting hole energy is 
-2y/3t. 



G 



Let us now come back to the original problem with the full Hamiltonian TCq +H.1 + TC2, Eq. (Q) and discuss how such 
an "improvement" of the Green's function proposed in Eq.(|l5|) can be done formally in the diagrammatic language. 

The fact that magnons in the string have lower energy than free excitations is the result of magnon-magnon binding. 
The magnon-magnon interaction (last term oiH.2, Eq. (|5|)) is attractive and leads to a bound state. A regular approach 
to the bound state problem is based on the Bethe-Salpeter equation shown in Fig. m 






+ 



FIG. 4. The Bethe-Salpeter equation for the two-magnon problem. Dots represent bare vertices T° b , and squares represent 
full vertices T a b- 

The equivalent integral equation is 

r ab (n, qi - q 3 ) = r« 6 ( qi - q3 ) + E r "fc(qi - Pi)r oi (fi, Pl - q 3 ) E z?»£>°(ft - w ) (16) 

Pl u» 

r o6 (n,pi - q 3 ) 



-4J7 qi -q 3 +^(-4j 7qi _ pl ; 



pi 



n - 4 j + zo 



It is easy to see that the momentum and frequency dependence of the vertex can be separated. The suggestion 
r a b(^,qi - q3) = -4J7 qi _ q3 f afc (0) transforms Eq. (|l|) into 



f «&(fi) = 1 



j 



(17) 



Here we have used the relation J2 Pl 7qi-pi7pi-q3 = 7qi-qs/4- Thus the renormalized magnon-magnon vertex is 
given by 



r a6 (n,k-p) = -4J 7k _ p 



n - 4 j + io 
n - 3 j + io 



(18) 



The pole of this function, f2 = 3 J, corresponds to the bound state whose energy is lower than that of two free magnons 
by — J, in agreement with our expectations. 

Similarly the hole- two-magnon interaction (second term in H2, Eq. (|J)) lowers the energy of the magnons in the 
neighborhood of a hole. The hole-three-magnon interaction (first term in 7i 2 ) originates from the projection operator 
and projects out the unphysical states with the hole and magnon at the same site created by the hopping term in TIq. 
A consistent consideration of these interactions yields the following Dyson equation for the single-hole propagator: 

G(e) = [ e -S( e )]- 1 , 

E(e) = E G(e - u)D\w) ]T r k _ q r k , q ( e , w ) , (19) 

uj q 

where the renormalized hole-magnon vertex is 

r k , q ( e , W ) = r°_ q + r« ( e ) + r< 2 q ( e , w ) + r k 3 q ( e ) + r£ q ( e ) , 

C(e)= E r°(q-c U )G( e - Wl )i3 (a; 1 )r kiqi (e,c I ; 1 ) , 

r k 2 q( £ - w ) = E E r k-q- qi ^ (^ - w)G(e - 0)r afc (tt, k - q 2 )^ ( W2 )^ (^ - ^ 2 )G(e - wj) 

qi,q2 W2,f2 

x r k . q+qi (e-w 2 ,r2-W2)r k ,q 2 (e,W2) > 

r k 3 q ( £ )= E E ^(k-q-q^ZA^ZA^Gte-Wi-c^ 
qi,q2 wi,W2 



(20) 
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r l 4 q( £ )= £ £ rO(k-q-Q)D ( Wl ) J D°(fi-c I ; 1 )G(e-r!)r Qb (r!,q 1 -q 2 )Z? ( W2 )C (r!- W2 ) 
qi,q2,Q wi,w2,n 

x r k , qi +Q(e - wa,n - Lo 2 )G(e - w 2 )r k , q2 (e, w 2 ) . 
A diagrammatic representation of Eqs. ([jl],|2C]) is shown in Fig. |[ 



+ 




FIG. 5. The Dyson equation for the hole Green's function. Bold lines are the dressed Green's functions, circles are the 
renormalized hole-magnon vertices, squares are the renormalized magnon-magnon vertices (Fig. ^| and Eq. dig)). 

Note that the lowest order correction to all vertices are zero, and higher order corrections are neglected. We have 
already discussed the absence of one-loop correction to the hole-magnon vertex. It is less evident in cases of the 
hole-two- and hole-three-magnon vertices. However, if one recalls the /, g and a, b indices for the hole and magnon 
lines in Fig. ^[ it becomes evident that the emission of an extra magnon prior to the hole-two-magnon interaction 
alters the type of the hole and therefore prevents it from absorbing and reemitting the original magnon. The same is 
true for the hole-three-magnon vertex. It is interesting to note that the set of diagrams in Fig. || is still in the SCBA, 
i.e., the crossing diagrams are absent. 

Again the momentum and frequency dependence of T can be separated. The almost evident substitution Fk, q (e, u>) = 
4t7k- q r(e,w) simplifies Eqs. ([l9 20) which, after integrating over the internal frequencies and momenta, becomes 



E(e) = 4t 2 G(e-2J)r(e,2J) , 
f(e,u) = 1 - G(e - 2J)f(e, 2 J) 



4r 



j 



UJ — J 



G(e - 3 J)f(e - 2 J, J) - G(e — 2 J — w)f(e - 2 J, u) 



(21) 



Each term in the square bracket comes from the corresponding irreducible diagram in Fig. || After combining similar 
terms T(e, 2 J) becomes 



+ t 2 G(e - 4 J)f(e - 2 J, 2 J) + t 2 G(e - 3 J)f (e - 2 J, J) - G(e - 4 J)f (e - 2 J, 2 J) 



f(e,2J) =G(e-2J)- 1 
= G(e-2J)- 1 



G(e - 2 J)" 1 + J/2 - 3< 2 G(e - 3 J)f (e - 2 J, J) + 4t 2 G(e - 4J)f (e - 2 J, 2 J) 



(22) 



e - 3 J/2 - 3t 2 G(e - 3 J)f(e - 2 J, J) 



G Q (e-3J/2) 



G(e-2J) ' 

where G a (e) has a continued fraction form. Summarizing Eqs. (|l9|) , (|2l|) , (p2[) one finally obtains 

G(e) = [e-4i 2 G (e-3J/2)]- 1 , 

ar, /, tm-1 _ 1 J-e/j(2y/3t/J) 



(23) 



G a (e) = [e-3i 2 G a (e- J)]" 



V3t J_ e /j-i(2V3i/J) ' 



which is exactly the same as the Green's function suggested in Eq. (|l5|). 

The rest of this Section is devoted to the comparison of the results from Eq. ( p3| ) with our numerical ED results 
on a 32-site lattice. Fig. ^| shows the ED results together with analytical results from the SCBA Eq. and from 
Eq. (H). 
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FIG. 6. The single-hole ground state energy v.s. t/J. The zero level of energy is taken to be the energy of the hole at t — 0. 

We recall that contributions from the Trugman's paths are left out from Eqs. (^3|) and (|^). These paths give rise to a 
small dispersion of the hole, SE^ = cos{k x ) cos(k y ) + (cos(2k x ) + cos(2k x )) + . . ., with t e ^ <C t\ . This form 
of dispersion implies that the correction to the energy at the point (tt/2,tt/2) due to the Trugman paths is almost 
zero. Therefore, it is natural to compare the ED data at this k-point to the analytical results. From Fig. ^ one can 
see a beautiful agreement of the numerical data at k = (tt/2, 7t/2) with the results of Eq. ( pq ) in the whole range of 
t/J. For comparison purpose, Fig. ^| also shows the numerical data at k = (0,0). We also remark that our 32-sites 
ED-data at k = (0, 0) agree with the 50-sites results of the modified Lanczos study by Riera and Dagotto (not 
shown in Fig. ^) up to the fourth digit at t/J = 5 and up to the seventh digit at t/J = 1.25. 

It is already clear from Fig. ^ that our present approach works better than the SCBA. To further illustrate this 
point, we plot the quasiparticle residue in Figure |7]. Three sets of data are shown. They are the results from our 
present analytical approach, the SCBA, and our numerical ED study. The ED results are calculated using the same 
method as in Ref. pU. One can see that such a comparison unambiguously favor the present approach. 



X 


- 32-site ED, k=(0,0) 


• 


- 32-site ED, k=(7i/2,3t/2) 




- - analytical result 




-SCBA 





FIG. 7. Quasiparticle residue for the lowest pole v.s. t/J. Dashed curve is the result of the SCBA. Solid curve is the result 
of the present approach. Dots are the numerical data for the (7r/2,7r/2) point. 
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V. TWO HOLES 



In this section we consider the problem of hole pairing in the t-J z model using the diagrammatic formalism in- 
troduced in previous sections. Our goal is to study in detail the bound state of two holes. The role of different 
interactions and higher order corrections will be analyzed. Bound states with s-, p-, and c£-symmetry and zero total 
momentum will be investigated. Throughout this paper we consider the S z = bound states only. Consequently we 
are interested in the interactions of the particles originating from different sublattices (/ and g fermions). 

If an exact two-particle Green's function has a pole in the scattering channel at fixed total energy of the particles, 
then there exists a bound state whose energy is the total energy of the particles. Thus, to analyze the bound state 
problem of two holes one has to solve the Bethe-Salpeter equation, whose most general form is given by 

r^ k , p 0W) =r£ kip (n, e,e')+ £ r^ kiPi (O,e,e 1 )Gp /2+Pl (O/2 + e 1 )G P/2 _ pi (n/2-e 1 )f^ i!p (Q, ei ,e0 , (24) 

Pl.ei 

where 

P = k!+k 2 =k 3 +k 4 = k 1 ,+k 2 , , k=(k!-k 2 )/2, p=(k 3 -k 4 )/2, Pl = (k r - k a ,)/2 , 

n = ei + e 2 = es + e4 = ei'+ea' , e = (ei - e 2 )/2 , e' = (e 3 - e 4 )/2 , ei = (e v - e v )/2 , (25) 

and the indices are given in accordance with Fig. [§. 

1 3 1 3 1 V 3 




FIG. 8. The Bethe-Salpeter equation for the two-hole problem. Lines represent renormalized Green's functions, dashed 
circles represent the bare "compact" vertex function F^ s , and black circles represent the exact vertex T . 

The bare vertex function V? 9 includes all diagrams which cannot be reduced to the second term of Eq. (Q), i.e., 
those which cannot be cut by a vertical line between the ends 1, 2 and 3, 4 into two parts joined only by two hole 
lines. The set of all such diagrams is called a "compact" vertex. 

At the pole, when f2 = E (E is the energy of the bound state), the left-hand side and the second term in the 
right-hand side of Eq. ( |2~4| ) are singular. Therefore T^ 9 ^> T^ 9 , and the first term in the right-hand side can be 
neglected. Hence the integral equation ( p4] ) becomes homogeneous in T^ 9 . The variables p and e' become parameters 
and are not defined by the equation itself. Omitting them and the indices /, g yields 

f P ,k(£,e) = J] r P)k)Pl (^,e J ei)Gp /a+Pl (B/2 + e 1 )Gp /a _ pl (B/2-ei)fp, Pi (£?,ei) . (26) 
pi>«i 

Further simplifications of this general formula Eq. ( pff ) are possible based on the specifics of the problem being 
considered. Rather general and frequently used simplification can be performed, for example, when the bare interaction 
r is a potential-like term, i.e. independent of the frequency. In this case T becomes e-independent and plays only 
an auxiliary role. Furthermore, it is more convenient to use the bound state wave function ipP.k(E) — GGT. By 
changing the variables and integrating both sides of Eq. (|26j) over e we obtain the Schrodinger equation for the bound 
state in momentum representation 

tfp,k(£) - ( Y, G P/2+k(£/2 + e)G P/2 _ k (i;/2 - e) J ]T rp, k , pi V>P, Pl (E) . (27) 

\ e / Pi 

If the Green's functions are the bare ones, the sum over e simply gives 1/[E — e(P/2 + k) — e(P/2 — k)], where e(k) 
is the single-particle energy. 

In the t-J z model the following characteristics allow one to simplify Eq. (|2^) considerably and obtain most of the 
results for the bound state problem in a transparent analytical form: 

(i) The single-hole Green's function is mostly k-independent, meaning that in a wide range of tj J ratio the k-dependent 
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contribution to the self-energy is insignificant or small. 

(ii) The lowest order contributions to the compact vertex rpk, P i {E, e, ei) have a simple kinematic structure and, 
together with (i), allow one to separate the momentum and frequency dependence of the exact vertex T. Its k- 
dependent part can be classified in terms of s-, p- and d-wave harmonics and can be integrated out easily. 

(iii) All one-loop and lowest order crossing corrections are exactly zero because of spin conservation in the holc- 
magnon interaction. Non-zero corrections are of higher order and are expected to be small. Therefore, the lowest 
order diagrams in T can be left unrenormalized. 

(iv) As we have shown in Sec. Ill the single-hole spectral function consists of a set of well separated quasiparticle 
peaks. Therefore, the Green's function of the hole can always be written in the form 



zO ' 



(28) 



where a u and e„ are the residue and the energy of the pole v, respectively. Then, if the characteristic binding energy 
is much smaller than the separation between the quasiparticle peaks (A = E — 2e a <C Se), i.e. one half of the bound 
state energy is close to the lowest peak (E/2 ps eo)> one can safely neglect the contribution of the higher poles and 
use the lowest pole approximation for the hole Green's function: 



G(e) 



«o 



e + iO 



(29) 



As a matter of fact, in the limiting case t J the separation of the levels is Se ~ J{tj 'J) 1 / 3 whereas the binding 
energy, as we will show, is E ~ —J(J/t). Thus, E/Se ~ (J/t) 4 / 3 <C 1 and the lowest pole approximation should work 
well. In the opposite limit i« J the binding energy is E ~ — J/2 while the separation of the quasiparticle poles is 
Se ~ 2J. Besides a factor of the order of E/5e ~ 1/4, the contribution of the higher poles is strongly reduced because 
their quasiparticle residues are negligible (~ i 2 /J 2 ), thus justifying the lowest pole approximation in this case. One 
can then argue that the same is true for all tj J. 

With these ideas in mind let us look at the compact vertex V. 




t t 



P/2-k P/2-p 




+ 



+ 




= r a +r b +r c 



(c) 



FIG. 9. The lowest order diagrams contributing to the "compact" vertex function F: (a) is due to nearest-neighbor hole-hole 
attraction, (b) are the one- magnon- exchange diagrams, (c) are the two-magnon exchange diagrams, (a), (b), and (c) are of 
first, second, and third order in Hint, respectively. Since the emission of magnon changes the spin of the hole, the positions of 
the ends of the lines in the diagrams in (b) are interchanged. P/2 ± k = (P/2 ± k, E/2 ± e), P/2 ±p= (P/2 ± p, E/2 ± ei). 

The lowest order contributions to the hole-hole compact vertex T are given in Fig. [| T a , r fc , and r c are of first, 
second, and third order in Hint =H\ +H2 (Eq. (^)), respectively. They can also be classified in terms of powers of 
J and t: T a - J, r fc - t 2 /J, and T c - t 2 / J . 

If we consider T b only, the full vertex T in Fig. ^ is given by the sum of the ladder diagrams only. To some extent, 
this approximation is analogous to the SCBA in the single-hole problem, and possesses an important property of the 
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latter: all crossing diagrams and one-loop corrections are exactly zero. More importantly, such corrections to all the 
lowest order diagrams (T a , T b , and r c ) are also zero because of the same spin-conservation rule. For the same reason, 
higher order corrections to Fig. ^ are strongly suppressed, so that the first dressing of T a , T b , and T c by magnons 
appears in the fifth, sixth, and seventh order, respectively. The corresponding set of diagrams is given in Fig. |l0|. 




FIG. 10. Lowest order corrections to the "compact" vertex function T (Fig. |9|). Circles represent T — T a + T b + T c , wavy 
lines are magnons. 

Note that because all these corrections are of similar origin to the Trugman's processes for the self-energy, one can 
expect that they play only a minor role in the effective interaction. 

However, there are other sources of correction to T coming from the hole-two-magnon, hole-three-magnon, and 
magnon-magnon interactions. Most of them result in the renormalization of the hole-magnon vertices similar to those 
in Fig. H While a formal account of such corrections is rather difficult, one can anticipate the result of such a 
renormalization using the real-space consideration of the interactions. As in the single-hole case, the dressing of the 
vertex accounts for two facts. First, the energy of a magnon in the neighborhood of a hole is — J/2 per hole-magnon 
link lower than the energy of a free excitation. Second, the coordination number for hole hopping in a state with a 
string of magnons is z — 1 instead of z . The first effect renormalizes the energy of the magnon lines in T b and F c , 
while the second effect renormalizes the inner hole line in T c . 

In what follows we will consider the compact vertex T as given by r a + T h + T c in Fig. || and will include the 
renormalization later. 

We first consider the t — limit when the solution of the BS equation (|26| ) is trivial. In this limit only T a — 
— 2J7k 1 -k 3 survives and the hole Green's function is simply 1/uj. Thus Eq. (pq) becomes: 

f P , t(E ,e) = -2JX>- P , g _^g^L — . (30, 

Evidently, the dependence of F on P is arbitrary and F(e) is a constant. Integrating over t\ one finds: 

2 7 

f k CE) = -— £>_ pi f pi (£), (31) 
pi 

where F plays the role of the wave function. Since the k-structure of the vertex is very simple, one can clas- 
sify all solutions of Eq. ( |3l| ) in terms of the nearest-neighbor s, p, and d waves: ip^. — cos(k x ) + cos(fc !/ ), 

(fk = V2 sm(k x ); -\/2 sin(fc y ), ^ = cos(fc.r) — cos(fc y ). Using the property: ^ p 7k-p¥'p P ' < ^ = j^k \ one can 
see that these solutions are degenerate with energy A s ' p,d — E s ' p,d — — J/2. This is expected if one remembers that 
while an individual hole breaks four AF bonds, two nearest-neighbor holes save one bond. 

A finite hopping constant t has three effects: (i) the hole Green's function becomes renormalized, (ii) r fc and r c 
start to contribute to the interaction, (iii) the hole self-energy also acquire a small k-dependent contribution. In order 



to see their roles, let us incorporate these changes into the BS equation (26) one by one. 

Consider the case when the hole Green's function is renormalized as described in Sec. Ill and is k-independcnt, 
but the interaction still comes from the T a term only. Since T a is P and e independent, the BS equation is similar to 
Eq. @ 

f k ( J B)=-2J^ 7 k_p 1 f pi (^)^G( J E;/2 + ei )G'(S/2-e 1 ) . (32) 

pi ei 

Since the k-structure of the core is unchanged, all three harmonics remain degenerate. Hence 
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1 = -^E G ( i? / 2 + e i) G ( £; / 2 - e i) 



(33) 



In the lowest-pole approximation G is given by Eq. ( p9| ) so the binding energy relative to the bottom of the hole band 
is: 



A = E- 2e = -Jal/2 



(34) 



In the limit t < J, A » -J(l - 16t 2 /9Jf)/2. In the limit i > J, a ~ J/i so A J(J/t) 2 . The binding energy 

(in units of J) v.s. f/J is shown in Fig. O. This energy is always negative because we assume that the hole forms a 
completely localized excitation. In this case a bound state is always formed as long as the effective interaction remains 
attractive. 

Now we consider the magnon exchange interaction T b which is given by: 



rLp(^i) 



7P/2-p7P/2-k 



7P/2+p7P/2+k 



ei-2J + iO -e-e 1 -2J + iO 



(35) 



Since T b depends on the total momentum of the pair P, the solutions of the BS equation also depends on P. The 
energy of the bound state at different P provides a dispersion law for the pair. In general, the lowest energy bound 
state is always at zero or some finite momenta corresponding to high symmetry points in the k-space. In this work 
we do not consider finite-P bound states because their physics is unimportant for the purposes of this paper. 
The BS equation @ at P = with interaction T a + T b reads: 



fo,k(-E, c) = -2jY,G(E/2 + et)G(E/2 - d) 



(36) 



8i_ 2 
1 



■7 P 7k 



ei - 2 J + iO 



ei - 2J + iO 



fo.p^ei) • 



Evidently, the k and e dependencies of f can be separated and one can see that the p- and d-wave solutions are 
not affected by the magnon-exchange interaction because J p fp^ = 0. This observation can be expressed in a 
different form: in the zero total momentum p- or d-wave state the amplitudes of the spin-flip exchange from different 
lobes of the wave function cancel each other. 

For the s-wave state, substituting fk(e) = tp^T(e) and integrating over k give: 



Me) 



J 



Y,G{E/2 + e 1 )G{E /2-eAl 



T 



i 



i 



e + e 1 -2J + iO -e-e 1 -2J + iO 



r B (ei) 



(37) 



In the so called static approximation, when the retardation due to the magnon propagator is neglected, T b is e- 
independent and Eq. (p7|) becomes 



(38) 



As in the previous case, A s = E s — 2eo- The problem with retardation can be solved using a method introduced in 
Ref. |24|j . One can assume that f b(e) has no singularities |^5| and therefore the integral over t\ in Eq. j37| ) is defined 
by the poles of G and the magnon propagator only. In the lowest-pole approximation, by using the obvious parity 
P(— e) = r(e) and by choosing the proper contour of integration one obtains: 



Me) 



2 A s 



1 



8f 2 



1 



1 



e + A s /2-2J + iO -e - A s /2 - 2 J + iO 



T E (A s /2) 



Multiplying both sides by G(E/2 + e)G(E/2 — e) and integrating over e finally yields: 



l = - J -f<l 
2 A s 



m 2 



J(A S - 2 J) 



(39) 



(40) 



This expression coincides with that given by the Wigner-Brillouin perturbation theory. The resulting binding energies 
(from T a + T b ) v.s. t/J are shown in Fig. |l| 
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FIG. 11. Binding energies (in units of J) of the s, p, and d bound states. When only the contact interaction T a is taken 
into account, all three bound states are degenerate (Eq. (p^)) and are shown by the lower curve. When the contact and 
magnon-exchange interactions F a + V b are considered, the p and d energies are unchanged (lower curve) while the energy of 
the s-state is pushed up (Eq. (ftOI)) (upper curve). 



Thus the s-wave binding energy vanishes at (t/J)\ c = l/(2v / 2) — 0.35. Since the spin-polarons are assumed to be 
localized the disappearance of the bound state is not due to the competition between the delocalizing kinetic energy 
and an attraction, but rather due to the competition between different types of interaction. One can see that the 
magnon exchange actually leads to an effective repulsion between the holes. In a naive point of view the holes can 
be mutually attracted by following each other in a caterpillar-like process, because in this way the disturbance on the 
AF background is only virtual. However, this scenario does not take into account the fact that after an elementary 
hopping of a hole pair, the order of the holes is different from the initial one, i.e., they are exchanged [ p6[ . This in 
turn makes the corresponding energy correction positive, i.e. it leads to repulsion. 

It is interesting to note that while the hole-magnon interaction plays a major role in the hole dressing, it has no 
direct effect on the p and d bound states and is destructive for the s-wave state. 

Let us now add the third part of the hole-hole interaction T c and see if the combination of the hole-magnon 
interaction with the two-magnon vertex causes an additional attraction. The diagrams in Fig. |j](c) are equivalent to: 



r£ )P (25, e, ei) = -2 J(4i) 2 7k _ p £ 7 £_ q £ —± 



u) - 2 J + iO 



G{e + E/2- u) 
uj - e + e x - 2 J + iO 



G(-e + E/2-u) 
u + e-e 1 -2J + iO 



(41) 



Integration over the inner momentum q and frequency u> yields: 



G(e + E/2 - 2J) - G(-e + E/2 - 2 J) - (e -> ei) 



(42) 



One can see that the kinematic structure of this vertex is identical to T a , i.e. this term is effectively a contact 
nearest-neighbor interaction. This is evident from the real-space picture of the process, which is shown in Fig. n3. 



o o 
oo-^ox^ox-^oo 

FIG. 12. Real-space image of the process in Fig. Blc) 
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Initially the holes are at nearest-neighbor sites. One of them jumps to a neighboring site creating a spin-flip at its site 
of origin. Then the second hole virtually absorbs and reemits this magnon. Finally the first hole jumps back and the 
original state is restored. Naturally, such an interaction depends only on 7k- P and does not depend on P. One can 
also guess the character of the resulting interaction. Being in the neighborhood of another hole, the spin excitation 
generated by the hopping of the first hole lowers its energy. Therefore the interaction between the holes via such a 
process must also be attractive. 

In the lowest-pole approximation r c Eq. (E2J) is given by: 



k, P 



(E, 



-8Ji 2 a o7k _ 



1 



(e + A/2 - 2 J + iO)(ei + A/2 - 2 J + iO) 



(43) 



One can consider the static limit of the problem by eliminating the frequency dependence of the propagators in 
Eq. This gives 



k, P 



-4— a 7k- 



(44) 



In the t ^> J limit this expression yields r£ ~ — 4i7k- p = — (2i/J)r£ p ^> r£ . Thus, in this limit the effective 
interaction constant of the contact interaction is renormalized from J to t. For the p- and <i-waves it gives a bound 



state energy A p ( d ) 



-tan 



-J(J/t), which is a factor of (t/J) 1 larger than the result obtained with T a only. 



This result, while being quite simple, is somewhat unexpected and, to our knowledge, has never been discussed in the 
literature. 

A rigorous account for the retardation yields modified equations for the bound state energies, leaving the qualitative 
picture drawn in the static limit essentially unchanged: 



1 = - J - 



2 ApW 



1 



8t 2 a 



1 = _Zfo 

2 A s 



(Ap(«Q - 2 J) 2 
8t 2 a 



p— , d— wave, 



(45) 



m 2 



(A S -2J) 2 J(A*-2J) 



s— wave 



Thus the energies of the p- and d-wave states are always negative. In the s-wave bound state, T c cannot overcome 
the magnon-exchange repulsion r fc . It also has little effect on {t/ J)\ s c because in the small t region where A s goes to 
zero, r c is roughly 4 times smaller than r fc . 

We now consider the effect of the renormalization. As we noted before we have enough reasons to rule out the 
possibility that the Trugman-like diagrams Fig. [Io| can change the interaction constants significantly. The anticipated 
smallness of these diagrams in the t > J region is due to the geometrical factor, which is of the order of l/z 3 = 1/64. 
Most of these corrections are relatively easy to evaluate using the same approximations as above. The results of these 
calculations show that those terms with the same kinematic structure as in r a , T b , and T c are indeed suppressed by 
this small factor. They also show that among the corrections there are terms with a more complicated kinematic 
structure than the original T. Since these terms are not suppressed by the small factor, they might enhance the 
effective interaction. Some of these terms correspond to the contact interaction between the holes in states where 
both of them have "tails" of strings. Evidently, this leads to interaction of polarons at sites farther than nearest 
neighbors. For the bound state wave function it would mean that harmonics higher than the nearest-neighbor one 
should be considered. Nevertheless, we do not expect these terms to make any qualitative change and in the following 
we will omit all corrections from Fig. ^ to avoid complications. 

At the same time, diagrams of the non-SCBA type (Fig. ||) are straightforward to include using the real-space 
language. We have to take into account the lower energy of the magnons in the intermediate state in the diagrams 
Fig. ||(b,c). This simply changes c^o = 2 J in the denominators of Eqs. (flq ) to ljq = J. Renormalization of the hole 
energy in the intermediate state in r c produces much weaker effect. This renormalization finally yields: 



J 

A s _ 

T 



< 1 



8t 2 a 



■fa 



St 2 a 



Jf 



(46) 



16t 2 



(A s -J) 2 J(A S -J) 



Solutions of these equations are plotted in Fig. 
changed significantly from those in Fig. [□]. 



O. One can see that the energies of the p and d bound states are 
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FIG. 13. Binding energies of the s, p, and d bound states when all three interactions shown in Fig. [] are taken into account 
and the energy of the magnons in the intermediate state is renormalized. 



Thus, the hierarchy of terms in the effective hole- hole interaction emerged from our study is as follows: 

(i) The contact "sharing common link" attraction yields bound states with s-, p-, d-symmetry with energy ~ — J/2 
at small t/J; 

(ii) The magnon-exchange interaction has no effect on the P = p- and d-wave states because of symmetry reasons, 
and is strongly repulsive for the s-wave state; 

(iii) The hole-magnon attraction, which leads to an attraction between a hole and a string of another hole, is the 
strongest interaction in the present problem in the physical limit (t > J). It leads to strong p and d bound states 
with energy ~ —J{J/t). 

Now, after all essential interactions are taken into account and carefully analyzed, the last effect which has to 
be incorporated into the BS equation before we can make comparison with numerical results is the k-dependent 
contributions to the hole Green's function. As noted in Sec. Ill such contributions are small, but lead to a coherent hole 
band. Therefore the ED data show some dispersion in the hole energy. For instance, at t/J = 5 the energy difference 
between the k = (0,0) and k = (vr/2,7r/2) points (almost exactly half of the bandwidth) is W/2 ~ 0.078t = 0.39 J. 
Although this difference is small compared to the absolute value of the lowest pole energy, eo ~ — 2.8t, or the separation 
between the energy levels, Se ~ t, it is as large as the hole-hole binding energy calculated above, A ~ —0.3 J. Therefore, 
one would expect this dispersion to have strong influence on the pairing of holes. One can also foresee that such an 
influence is likely to be destructive, because now the pairing is not only an interplay between repulsions and attractions, 
but also the question of how much kinetic energy is lost if the holes are in a bound state. 

For our purposes we will need the hole Green's function in the lowest-pole approximation, which now reads: 



G(e) 



«o 



eo 



Wt^/2 + iO 



(47) 



where £q — e — W/2, e k = (1 — cos^) cos(fc y )) > 0, W > 0, and we have neglected the k-dependence of a . The 
role of the Trugman's loops larger than an elementary square is negligible, so we consider the dispersion given by the 
next-nearest-neighbor hopping only. The resulting band is isotropic and has a minimum at the center of the BZ. Now 
in the BS equation fl2q ) it is easier to integrate over t\ first. Neglecting the effect of the dispersion on the retardation 
and after some algebra one obtains: 



= -2Jal^ 



9p 



^ A 
p 



lFe r 



7k- 



1 



8t 2 a Q 
(AT~7p 



m 2 



■7k7 P 



J(A - J) 



(48) 



where ip^ is a nearest-neighbor s-, p-, or d-wave. As before the magnon-exchange interaction (the second term in the 
curly bracket) is orthogonal to the p- and d-wave states. Eq. ( |48| ) is an almost evident generalization of Eqs. ( |4^ ) 
to the case of finite hole dispersion. The only new parameter introduced is the hole bandwidth W. Our numerical 
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ED data for W/2t as a function of t/J calculated on a 32-site lattice are given in Fig. We do not calculate this 
quantity analytically. Instead we treat it as a parameter and use the numerical values in Fig hj. 



0.10 



0.08 



0.06 



0.04 



0.02 



0.00 




FIG. 14. Half of the hole bandwidth in units of t as the function of t/J. 

Two statements can be made about the solutions of Eq. (f48|). Now all three symmetries have thresholds (t / J)\^.' p ' d 
above which negative energy solution of the corresponding symmetry cease to exist. Also, the p- and d- waves are split 
with the p wave state having lower energy. 

The first statement is evident if one consider the threshold condition (A = 0) for the p- or d-wave state in Eq. (El), 
assuming t ^> J so that ao — J/t: 



1 ^ --^r 



2 i 2 ^ -We* 



8t 2 Jl ^ 47 J 



^d 



(1 — cos{k x ) cos(ky)) 



(49) 



The sum over p gives a number of the order of unity. Thus, there is a value of W ((W/t) c ~ 4 J 2 /t 2 <C 1) at which the 
equality condition is satisfied. If W grows with t, as it in fact does, then above this point Eq. p§) has no solution. 
Using the dependence of W on t/J (Fig. |lj), the critical values of t/ J\ p,d can be extracted from the same condition 
Eq- (II). ' ' 

The second statement is a consequence of the isotropic form of the dispersion law. The sum over p in Eq. (|49| ) can 
be rewritten as: 



A p,< 



,nP,d 



ip?p ,d [l + cos(fc x ) cos(k y ) + (cos(k x ) cos(A; y )) 5 



■] = 1 + Ci + C 2 



(50) 



In the case of d wave the series is alternating with C\ — —1/2 and C2 = 3/8. In the case of p wave C n = for odd n, 
and C n > for even n. Numerical integration gives AP = 1.27, ^ d = 0.73, and therefore d wave solution disappears 
at a lower value of t/J. 



Numerical solutions of Eq. (48) for all three symmetries together with results of our 32-site ED calculations are 
shown in Fig. [l^. For the purpose of comparison we also show the results obtained by a modified Lanczos method 
on a 50-site cluster in Ref. 1231. 
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FIG. 15. Binding energies of the s-, p-, and d-wave bound states. Bold lines are the solutions of Eq. (w8h- Triangles, 
diamonds, and solid circles are the 32-site ED results for the s-, d- and p-wave, respectively. Lines connecting these points are 
guides for the eye. Empty circles are the results of modified Lanczos study on a 50-site cluster from Ref. |23|. 



One can see a very good quantitative agreement between our analytical and numerical results for the s-wave state. 
For the d- and p-states, they have the same qualitative behaviors. It is clear that our Eq. (|48f), while involving 
different kinds of approximations, captures all essential effects of the pairing problem in the t-J z model in a very 
good quantitative level. It also describes the problem in a very easy and transparent way. Namely, all essential 
tendencies in the hole binding are described by this equation: (i) the s-wave is pushed up to the continuum by the 
magnon-exchange process, (ii) the p- and d-waves are much more stable, (iii) the Trugman's dispersion split these two 
levels and push them to the continuum at some larger critical values of t/J. 

The analytical results underestimate the pairing strength, especially in the region i/J ~ 1. We believe that this is 
due to the contribution from the higher poles which are neglected in our study. In this region a significant portion of 
the spectral weight is transfered to these poles while the separation between them is still of the order of 2 J. Since 
our lowest-pole approximation uses the ratio A/Se (~ 1/4 in this region) as the small parameter, one would expect 
the energy correction from the higher poles to be of the same order. In the region where the binding energy is close 
to zero we expect this approximation to be more accurate. However, for the p and d waves this is also the region 
of t/J > 2 where the omitted higher order diagrams with more complicated kinematic structure contribute to the 
interaction. We refer the smaller splitting of the p- and d-wave states found in our analytical results compared to 
those in our ED data to the effect of these corrections. The same higher order interactions cause a tiny splitting of the 
states at t/J < 1 which is shown in Fig. [lf| We remark that a t/J dependence of the binding energies very similar 
to our ED results has been found recently by series-expansion calculations |2^]. Thus we believe that the agreement 
between the results can be improved by an accurate consideration of the contribution of higher poles and higher order 
diagrams of various types. However, this is out of the scope of our paper. 

It is natural and instructive to consider the case where small transverse spin fluctuations J± = a J (with a <C 1) 
are included in the Hamiltonian. This helps one to see the effects of these fluctuations on the pairing problem and to 
infer the tendencies in the isotropic t-J model. As noted in Ref. |l7|] only two changes are necessary to account for 
these fluctuations up to first order in a. First, the holes acquire a dispersion law given by 

5e k = aA^l (51) 

where A is of the order oft 2 / J in the t <C J limit and of the order of J in the t ^> J limit jL7|. Such a band is strongly 
anisotropic, i.e. degenerate minima form a line along the magnetic BZ boundary. In the full t-J model (a — 1) this 
dispersion dominates over the Trugman's terms at all t/J. Second, the hole-magnon interaction becomes: 

r S = 4i[7k- q - ^«7k7q + 0(a 2 )} . (52) 
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Corrections to the magnon energy are of higher order in a: luq = 2 J + 0(a 2 ). As we shall see both changes lead to 
the splitting of the p- and d-w&ve levels, and both favor the d-wave as the ground state of the system. 

Consider Eq. ([l8|) with the dispersion from Eq. (^lj) and the hole-magnon exchange term coming from Eq. (S3). 
For the purpose of illustration it is sufficient to consider the interactions in the static limit: 



p p 



1 



8t 2 a 



16t 2 
J 2 



7k7 P - 2^ k +p(Tk + 7 P ) 



(53) 



Corrections to the bound state energies of the p and d wave up to first order in a are: 



-a, 



where the first term in the bracket comes from the dispersion and the second is from the magnon-exchange interaction. 
The kinetic energy pushes up the energy of both states, but for the p-wave this shift is three times larger. An important 
observation is that the magnon-exchange interaction is now attractive for the <i-wave and repulsive for the p-wave 
state! 

To show how subtle the selection between the p- and d-wave ground state is, let us consider a more "realistic" 
dispersion law: 

<5e k = ^A(cos(k x ) 2 + cos(k y ) 2 ) . (55) 

It is more realistic in the sense that it has minima at ±(±7r/2, 7r/2) points which are isotropic and it agrees with 
the experimentally measured profile in the cuprates better than the one in Eq.(^). It is now known that this 
difference between the t-J model dispersion law (Eq.©) and the one observed in reality (Eq. ( |5^)) is due to additional 
intrasublattice hoppings (t' ,t" etc.) in the real systems ||. With the dispersion law in Eq. (|55D corrections to the 
binding energies are given by: 

A rf = A^(W^) , (56) 
A* = AS + |(3A+^ 

Now the kinetic energy pushes up the <i-wave energy 5/3 times more than that of the p-wave and A p — A d = 
j[2t 2 aQ/J — A] is not necessarily positive. Therefore, p-wave can become a ground state again. 

Summarizing our observations in the small a limit, we conclude that the transverse fluctuations create an additional 
attractive interaction in the ci-symmetry state and that an anisotropic kinetic energy profile also favors <i-wave. Thus, 
we suggest that the d-wave is most likely to be the ground state of the system in the t-J model (a — 1). This conclusion 
is in agreement with many analytical and numerical studies. Nevertheless, recent 32-site ED study showed that at 
tj J > 3.3 the ground state of the two-hole t-J model is a p-wave state with total momentum P = (ir, it) p7|. (In our 
case momenta P = (0,0) and P = (7r,7r) are equivalent because of the existence of long-range order.) Furthermore, 
such a behavior is also suggested by recent series-expansion calculations [^8| . Due to the apparently strong finite-size 
effects on the ED results, it is difficult to conclude whether this behavior will survive in the thermodynamic limit. 
Nevertheless, these new data show that even in the pure t-J model the final word on the problem of the symmetry of 
the bound state is still to come. 

Besides effects arising from the t-J model there are also concerns about the robustness of the bound states in the 
real system. The characteristic bound-state energy found in the t-J model is of the order of — 0.1J in the physical 
region of tj J. This energy is far below the accuracy of the model itself. This means that in order to describe the 
real materials, the low-energy model must include at least additional hopping terms {t',t" etc.). It is worth noting 
again that because of these terms the realistic dispersion is isotropic, which, as we have just shown, can favor p-wave. 
More important conjecture can be made using the results of Ref. Q and new ED data Since the isotropic 

dispersion involves higher kinetic energy, the overall pairing tendency is weakened by these additional hopping terms. 
It can push all the bound states up to the continuum, thus destroying the foundation of the preformed pair and phase 
separation scenarios well before the physical region of parameters is reached. 
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VI. CONCLUSIONS 



In this paper we attempt to address the issue of how generic d-wave pairing in t- J-like models is and to clarify 
the problem of how this symmetry is selected among others. We have thoroughly examined the t-J z model by means 
of analytical diagrammatic study supplemented by our 32-site ED results. This formalism is tested in the one-hole 
problem and perfect agreement with ED results is found. The ground state of two holes in the physical region t/J>l 
is found to be a p-wave bound state. From a careful analysis of the interactions involved in the pairing problem 
we established that the pairing in the s-wave channel is highly unfavorable because of the repulsive character of the 
magnon- exchange interaction in this state. While this reason is rather generic, the problem of selecting between the 
p- and d-wave is much more subtle. In our analytical study these states remain degenerate until the hole dispersion 
is included. In the t-J z model such a dispersion comes from higher order hopping processes and is, in fact, very 
small. However, it has important consequences. It is responsible not only for the significant splitting between the p- 
and d-wave states, but also for the destruction of hole pairing at some critical (t/J)\P' d . The isotropic form of the 
dispersion favors p-wave as the ground state of the system. 

We have also considered the role of transverse fluctuations in the perturbative limit a <C 1. They lead to an 
anisotropic hole dispersion and an additional hole-hole interaction. We have shown that both features favor d-wave as 
the lowest energy state. However, in this case the selection between p- and <i-wave strongly relies on the anisotropic 
form of the hole dispersion, and such a dispersion does not agree with the experimentally measured dispersion in 
real systems. Moreover, since the characteristic energy of the bound states found in the t-J models is small, we are 
skeptical that these states can survive when additional kinetic energy terms are included in the model. Our skepticism 
is supported by earlier studies and recent numerical data |29|j3C| ]. 

Such conclusions for the pairing problem in t- J-like models drawn from our study create a framework in which 
it might be necessary to look for d-wave superconductivity unrelated to the Bose-condensation of preformed d-wave 
pairs. Furthermore, to make the phase separation scenario possible, one might need to invoke arguments other than 
an a priori existing collapsing tendency in the spin-hole system. 
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